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' We derive a simple closed form for the matrix elements of the quantum baker's map that 

, shows that the map is an approximate shift in a symbolic representation based on discrete phase 

0^ ' space. We use this result to give a formal proof that the quantum baker's map approaches a 

classical Bernoulli shift in the limit of a small effective Planck's constant. 

C : 1 Introduction 

O^' The quantum baker's map [0, ^] is a prototypical quantum map invented for the theoretical study of 

quantum chaos. During the last decade, its semiclassical properties have been studied extensively 



^ ! 1^) B 1> H 0) 1] ; it has been shown to display hypersensitivity to perturbation |T^, optical |jTl|] 

' and quantum computing |l^, |l^ realizations have been proposed, its long-time behavior has been 

investigated @, 0, it has been studied in a path-integral approach |^] and defined on a sphere |14|. 



The quantum baker's map is a quantized version of the classical baker's transformation ||T^, but 
there is no unique quantization procedure |16|. The original definition of the map |l|, |2| is based on 
Weyl's quantization [^] of the unit square. Essentially the same map has been derived by algebraic 
methods [18, 19 1 as vi^ell as by considering the transition from ray to wave optics [|Tl|]. Recently a 



whole class of quantum baker's maps has been defined |20| by exploiting formal similarities between 
the symbolic dynamics pl| ] for the classical map on the one hand and the dynamics of strings of 
quantum bits of the type considered in the theory of quantum computing on the other hand. This 
class of quantum baker's maps, which can also be derived from the semiquantum maps introduced 
in Ref. |^, is the subject of this paper. 
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^ r.schack@rhul.ac.uk 



1 



The classical baker's transformation [|T5| maps the unit square < q,p < 1 onto itself according 

to 

, , I (29,IP) , ifO<g<i, 

{q,p) 1 ^ < > 1 , ,\ , (1) 

[ (2g-l,i(p+l)j , if i <g< 1. 

This corresponds to compressing the unit square in the p direction and stretching it in the q 
direction, while preserving the area, then cutting it vertically, and finally stacking the right part 
on top of the left part — in analogy to the way a baker kneads dough. The classical baker's map, 
has a simple description in terms of its symbolic dynamics |2l|]. Each point {q,p) is represented by 
a symbolic string 



' S-2S-1S0.S1S2 ■■■ , (2) 



where Sfc = or 1, and 



g = ^.fc2-^ p = j;._fc2-^-i. (3) 

k=l k=0 

The action of the baker's map on a symbolic string s is then given by the shift map (or Bernoulli 
shift) U defined by Us = s', where = s^+i. This means that, at each time step, the entire string is 
shifted one place to the left while the dot remains fixed. Although the relation (|3|) between points 
{q,p) and symbolic strings is particular to the baker's transformation, the method of symbolic 
dynamics is very general and can be applied to a large class of chaotic maps pl| . 

Symbolic representations for the quantum baker's map have been introduced in Refs. [|, pi, |20] . 



These representations are all obtained by writing the quantum propagator B in a mixed form 
{i\B\i'), where {\i)} and are different bases. In this paper, we derive a simple closed form 

of the matrix elements with respect to a single basis. We show that all members of the class of 
quantum baker's maps defined in Ref. are approximate shifts in a symbolic representation 
based on discrete phase space. We use this result to give a formal proof that all members of this 
class of quantum baker's maps approach a classical Bernoulli shift in the limit of a small effective 
Planck's constant. 

The paper is organized as follows. In Sec. we give the necessary background and definitions. 
In Sec. ^, we state the results of the paper and discuss their significance. Finally, Sec. |^ contains 
the derivations and proofs. 



2 Background 

Most results of this paper are phrased in terms of finite binary strings. It will be convenient to 
adopt a slightly different and more flexible notation than the one used in Eq. (|2|). Here, a binary 
string 

def / 66+1 • • • (s < /) / 

where 6 £ {0) 1} is a bit, can have increasing {s < f) or decreasing (s > f) indices. We will use 
bold Greek and Latin letters to denote binary strings, e.g., 

a = 6:/ or a; = Xh-.t ■ (5) 
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The length of a string a wih be denoted by \a\; e.g., in the above example, \a\ = |/ — s| + 1. 
Concatenation of strings is defined in the usual way. Again considering the above example, ax 
is the string cxx = S^s---(,fXh---Xt- Any string a represents a natural number through its binary 
expansion 

|«| 

a = E2'"'"Sfc)' (6) 

k=l 

where a^^) denotes the A;-th bit of a, 1 < < |a|, such that 

a = 0(1)0(2). ..Q!(|q,|) . (7) 

Thus our notation does not distinguish between a binary string and the corresponding natural 
number. Similarly, two strings a and x can be combined to represent a rational number 

a.x 2l"l-'«W + E • (8) 

k=l k=l 

Quantum baker's maps are defined on the Z)-dimensional Hilbert space of the quantized unit 
square [17|. For consistency of units, we let the quantum scale on "phase space" be 2TTh = 1/D. 



Following Ref. Q, we choose half-integer eigenvalues qj = {j + ^)/D, j = 0, . . . , D — 1, and 
Pk = {k + ^) / D , k = 0, . . . ,D — 1, of the discrete "position" and "momentum" operators q and p, 
respectively, corresponding to antiperiodic boundary conditions. We further assume that D = 2^, 
which is the dimension of the Hilbert space of N qubits, i.e., N two-state systems. 

The D = 2^ dimensional Hilbert space modeling the unit square can be realized as the product 
space of N qubits in such a way that 

kj) = 16)55 16) ®---®iejv), (9) 

where j = "^^iLi 6 ^ {Oi l}i ™d where each qubit has basis states |0) and |1). It follows that, 

written in our string notation as binary numbers, j = ■ ■ - ^.N = (.i-.N and qj = O.S,iS,2 ■ ■ ■ Cn^ = 
O.Ci:Afl- We define the notation 

|.ei:^) = |.66...e^) = e^"/'k,), (10) 

which is closely analogous to Eq. (|2|), where the bits to the right of the dot specify the position 
variable; see Ref. |2^] for the reason for the phase shift e*"^/^. 

Momentum and position eigenstates are related through the quantum Fourier transform oper- 
ator F Q, i.e., F\qi:) = \pk)- Again in analogy to Eq. (|2|), we define the notation |6:Af-) = \Pk)i 
where pk = O.^Ar:il. 

By applying a partial quantum Fourier transform to the position eigenstates, one obtains 
the family of states 

\Cl:n.Cn+l:N) = \Cn+l) ® • • • ® |e7v)e'"(°-«-^^)® 

v/T72{|0) +exp[27ri(0.eil)]|l)}® 

Vl/2{\0) + exp[2^i(0.661)]|l)} «>•••«> 

^/l/2{\0) + exp[2^i(0.e„;il)]|l)} , (11) 
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where 1 < n < — 1. More precisely, the state \Ci:n-S,n+i:N) is obtained by applying the Fourier 
transform operator to the n rightmost bits of the position eigenstate |.^n+i:AfCn:i)- For given n, 
these states form an orthogonal basis. The state \£,i:n-(,n+i:N) is localized in both position and 
momentum: it is strictly localized within a position region of width l/2^~", centered at position 
Q = 0.^n+i:Afl; and it is crudely localized within a momentum region of width 1/2", centered at 
momentum p = 0.i^„:il. 

For each n, l<n<N — 2, a quantum baker's map can be defined by 

-B|^l:n-'^n+l:Af) '= \(.l:n+l-S,n+2:N) , (12) 

where the dot is shifted by one position. In phase-space language, the map B takes a state local- 
ized at {q,p) = (0.^„+i:Arl, 0.^n:il) to a State localized at {q',p') = (0.^n+2:Af 1, 0.^n+i:il), while it 
stretches the state by a factor of two in the q direction and squeezes it by a factor of two in the p 
direction. This analogy with the classical baker's map motivates calling the maps ( [T2|) "quantum 
baker's maps." For n = N — 1, the map is the original quantum baker's map as defined in Ref. Q, 
which in our notation becomes 

mi:N-UN) = \Cl:N.) , (13) 



and for n = 0, the map is 



BUi:N) = \^U2:N) ■ (14) 



Below we show that all the maps (12,1^,14) reduce to the classical baker's map in the limit h ^ 0. 



3 Results 



Equation (12) is a mixed representation of the quantum baker's map, using different bases on both 
sides of the equation. To go beyond the heuristic phase-space interpretation of the map given at 
the end of the last section, we need to express the matrix elements of B with respect to a single 
basis, i.e., we need to find 

^Ist/tO dcf I {■(,1:n\B\-S,1;n) if 71 = 



where e=^U and ^=^1 

.jy. A main result of this paper is the following simple formula, which 

will be proved in Sec. 0: 



^ '^) = ^(^-^-^2-^1 sin[vr(0.e°,,,l- 0.4,1)]' ^''^ 



where 1 < n < N — 2 and <^ is a phase factor given by 



c^(e?,d) = ^w-i)«--(-i)«"]- (17) 



For the case n = 0, one obtains 



C'^\e,e) = "-^^iCU - Cj:iV-i)e^^'«?-«-l , (18) 
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and for n = N — 1, 



N) 



(19) 



The coefficients C^^^{^^ ^ given in (16) are zero unless ^^_|_2 Ar = '^j^+i-at-d i-e-, unless the position 
bits Cre+i AT-i of the final state are obtained by shifting the corresponding position bits of the 
initial state. Furthermore, the sin term in the denominator ensures that C^^^ is strongly peaked 
for £,n+i-2 ~ '^n i' i-^-' if ^be momentum bits ^^.^ of the final state are obtained by shifting the 



corresponding momentum bits of the initial state. The formula ( fig) therefore establishes rigorously 
that the maps (|l^ ) are approximate shift maps, a result which had been obtained numerically in 
Ref. i. 

To formulate the question of the classical limit of the baker's map, we use the concept of coarse- 



graining in the spirit of the consistent (or decoherent) histories approach |22, 24 1. For this. 



we introduce projectors on subspaces corresponding to symbolic strings y of length I. We fix in 
advance an upper limit, /cmax) on the number of iterations, k, considered; this is necessary because 
in computing the classical limit of a chaotic map, the limit h = 2^^^^^^ /tt has to be taken 
before the limit /c ^ oo |25|. We will show that, for given / and k^^ax, it is always possible to choose 
h in such a way that the coarse-grained quantum dynamics is arbitrarily close to a shift of the 
string y. In contrast to the approach of Refs. p^, 19|, in taking the limit ^ ^ 0, we always remain 



in the finite-dimensional Hilbert space on which our maps are defined. 

As before, we are considering basis states of the form \£,i:n-Cn+i:N) ■ As we let N increase, the 

def 

number of position bits to the right of the dot, m = N — n, remains fixed. We define r = N — I 
as the number of bits ignored in the coarse graining. In the following, we always assume that 

k ^ fcmax < 

We are now in a position to introduce a family of projectors 



^ \xy^.y^g){xy^.y'^g\ if k < 



pr,k <l|f 



m 



\x\=r—k. \g\=k 



^ \xyg^.g^){xyg^.g^\ \ik>m, 



(20) 



|a3|=r-fc, |g2 

I —k — m 



where y^y^ = y and \y'^\ = m — k. By normalizing these projectors, we obtain a family of uniform 
density matrices. 



A classical shift acts on these states as 



Pk 1-^ pk+i 



Projecting a state pk' onto the shifted subspace Py'^ gives the characteristic delta distribution 



We will prove that 



Tr [P^'Vfc'] = hk' 



(21) 
(22) 

I 

(23) 
(24) 
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or, since k is bounded from above by fcmaxi and r = N — I, where I is fixed, 

Tr [P^''=^Vo(^^)'] = 1 - O(^) = 1 - 0{hlogh) . (25) 

Comparing Eqs. ( |23| ) and (25), one sees that the coarse-grained quantum evolution approaches the 
shift-map behavior to any required accuracy as h ^ 0. A measurement of the projectors Py can 
be interpreted as a measurement in which the r — k leftmost bits and the k rightmost bits of the 
symbolic string are not resolved. 

Equation ( p^ ) can be rewritten as 

2- Tr [P^''B'\xy\y^){xy\y\B^)''] = 1 - O(^) , (26) 

\x\=r 

which is a sum of 2^ terms bounded from above as 

Tr [P'^'''B''\xy\y^){xy\y\l3^)''] < 1 . (27) 



Here, y^y^ = y and = m as before. Equations ( p6|) and (p7|) can be both satisfied only if the 
condition 

Tr [P;^''B^\xy\y^){xy\y\B^)'^] = 1 - O(^) (28) 

holds for all x except for a fraction of order r/2'^~^, i.e., for all basis states \Ci:n-Cn+i:N) except for 
an exponentially small fraction. In other words, the property (^) holds for typical basis states. 

An interesting feature of the quantum baker's map is that there are atypical basis states for 
which Eq. ( p8[ ) does not hold. In section IV we give an example of an atypical state \x^^^^y^ .y"^) 
for which 

Tr [P;^'^B\x''yPy\y^){x''yPy\y^\B^] = ^^^^^ + 0(4'^") + 0(2"^) , (29) 



where G ~ 0.915965 is Catalan's constant [26|. For sufficiently large n — r and r, this expression 
is less than 0.872. This is an example where the quantum evolution in the limit — > differs 
substantially from the classical evolution, already after the first iteration of the map. If, however, 
the initial state is a mixture in which atypical states have an exponentially small weight, such as 
po, the correspondence principle is obeyed. 



4 Derivations and proofs 
4.1 First Iteration 

In this section we prove the formula ( |T^ ) for the matrix elements C^'^*(^'^, ^^). Equation (|^) for 
the case n = — 1 follows from almost identical arguments, and Eq. ( p!8|) for n = is essentially 
trivial. A direct calculation yields 

c'^\e,e) = 5(e°+2-e^+i)^(e+3-e^+2)---5(d-d-i) 

x,/T/2{6{eN)+HfN - l)exp[2^i(0.e?l)]} 
xexp{m(0.eO+i^il-0.4il)} 
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xl/2[l + exp{27ri(O.C5C?l - O.^Jl)}] x ••• 

xl/2[l + exp{2^i(0.e°+i:il - O.^il)}] . (30) 

Using the identity I + e^'f> = 2e"^l'^ cos((/>/2) and noticing that (5(^j^) + - 1) exp[27ri(0.^?l)] 
=exp[i7r^^(^]' + 1/2)], we have 

C^^\i\e) = yi72exp[ivre^(e? + l/2)]5(e°+2:7V-^n+l:^-l) 

xexp[i7r(0.e°+i;il-0.eiil)] 

n+1 

X jjcosKo.eLi-o.eLi^ii)] 

k=2 
n+1 

xnexpM04°il-0.eLi;il)] 

fc=2 

= yT72exp[i7re^(e? + l/2)mn+2:N - Cl+l:N-l) 
n n 

xe'^-{\[co^(t>k){\{e'^'^) , (31) 

k=l k=l 

where 

./-fc = vr(0.4Vil-0-4il)- (32) 

To simphfy Eq. (^Tj), we first consider the products of cosines and exponents separately and 
then combine them to formulate the final result for the first iteration of the baker's map. Note that 

24>u = 7:{{).ilil-Ul-i:il)+<ek+i-il) 

= 0,_i + ^(eo+,-ei), (33) 

so 

cos 4)k-i = cos[2</)fc + ■n{il - il^i)] 

= {-ifl-^l+i cos{2(f)k) , k<n . (34) 

From Eq. (32), we have 2(f>k = (mod 27r) and thus 2(j)k = 2"^^^'^'0„ (mod 27r), so the previous 

formula can be rewritten as 

cos(j)k = (-l)^^-+i"^"+2cos(2"-''(/.„) , k<n-l. (35) 

Using this formula the product of cosines can be expressed as 

n n— 1 



cos (pk = COS (pn Y\_ COS 4>k 

k=l 

n 

JJ cos[2-'=(2"0„)] , (36) 



k=l k=l 



k=l 



where a{£,k:n) '= X]s=A:^«- ^^^^ ^° check by induction that 



Tl 1 . r\fi 

nu sill Ax, , 

cos2'^x = — x/vrj, j = 0,±l,±2,... 37 
2" sm X 

k=o 
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In our case 

k=l 

Putting everything together, the product of cosines becomes 

n = (-1)^(51: J-(?L...) , (39) 



2" sin ( 

k=l 



where 4>n = "O-'^n i-'-)- Now we simphfy the product of exponents in (|3l|). Equation (^3|) 

imphes 

k 

2Vn + 5^ 2'=-^7r(ei+i_, - e°+2- J, ^ > 1 , (40) 



so 



n-1 



fc=l 
n— 1 n— ] 



k=l k=l 

n—1 n—1 k 



k=0 k=l s=l 

n-1 



= [<An(2" - 1) + ^7r{en+i.k - en+2-fe)](mod 27r) 

k=l 

= {<A„(2'^-l)+^[a(CL)-^(e3°n+i)]}(mod27r). (41) 
The product of exponents is thus given by 

fr e'^'^ = exp(i y <t>u) = (-l)-(?2. J-«L+i)!^Pi!^ . (42) 



Using (|39|) and (42) one can rewrite (31) as 

,u _ (-l)g^(g?+V^) 1 sin(2"0O 

<^ ^„^,i/2 ^\^n+2:N - C,n+l:N-l)—^^ exp(i2 0„) . (43) 

Further simphfication is possible due to the fact that 2"(;6„ = 27r(0.^2?]'l — 0.^|l) (mod 27r). The 
final result is 

^ ^^'^) = ^(^-^-^2"+isin[vr(0.eVil- 0-4.1)]' ^''^ 
where the phase factor <I> is given by 

ml eir) = \/2e^-«^«?+^/2) sin(2>,) exp(i2>„) 

= -L[i(_i)C^ _(_!)«?]. (45) 

This formula is an exact expression for the matrix elements ( [l5| ) of the quantum baker's map. 
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4.2 k-th iteration 

In this section we prove that for all k < k^aax^ 

jr,k ^k^^/nUk 



r 



Tr[P;''B'poiB^r] = l-0{^), (46) 



where the projectors Py and the density operators pj are defined in Eqs. ( pO[ ) and (pi]). The first 
step is to prove that 

Tr [P^'^'Bp,_iBt] = 1 _ o(^) . (47) 

By a direct calculation, we obtain 

TT[P^'''Bpk-iB^] = 2^-'- |2''+^sin[2-("-''-+'=)7r(0.al-0./31)]r2 

Ct|=r-fe + l 

^ E E - 4- - 1)-' ' (48) 

n=0 D=0 

where the inequality sin^ x < was used. Let L = r — k, and let Q(s) be the number of different 
pairs {u,v), < u < 2^+^, < v < 2^, for which u — 2v = s. It follows that 

E E (2n-4T;-l)-2 = ^ Q(s)(2.-1)-^ 

u=0 v=0 s=_2(2i-l) 

''^'' Q(g) + Q(i-^) 

- (2s -1)2 ■ ^^^^ 

s=l 

Using a simple counting argument based on the register principle |28], one can show that 

Q{s) + Q{1 -s) = 2^+1 - s + ^[1 - (-1)1 , (50) 



Q(s) + Q(l-s) 2^+1 -s 



from which one obtains 

L+l 

^ (2.-1)2 ^(2.-1)2 

4 4^ 



(2. - 1)2 ^ (2. 

^"^[4 - Oiir)] , (51) 



where we have used the relations 



E(2^ - 1)-^ = Y + 0(2-^) , E ^ = O(^) . (52) 



s=l s=l 



Combining Eqs. (^), (|4^ ) and (|5l|), we obtain Eq. (47) as required. We now rewrite Eq. (^) in 
the symmetric form 

Tr [pkBpk^.B^] = 2-^[l - 0(-^)] (53) 
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and introduce the distance measure between density matrices induced by the Eucledian norm 

= VTr(p-p')2. (54) 

This distance measure is unitarily invariant and obeys the triangle inequahty. We will now prove 
that (H) and (|5§) imply 

d{pi,,B^po[B^f) = 0{2'^-'-^) . (55) 
Using the cyclic property of the trace, we have 

d{pk, Bpk-iB^) = ^^^pI + Tr pl_^ - 21V {pkBpk-iB^) . (56) 
Since Tr pi = 2^ jl^^ = for any k, 

[d{pk,Bpk-iB^)f = 2-^+^ - 2Tr {pkBpk-iB^) , (57) 
which, together with Eq. (|53|), implies 

d{pk,Bpi,_iB^)=0{2-2-^V^) . (58) 



The case k = 1 of (|55| ) follows directly from ( |58D . Assuming that (|5^) is true for a given value of k 
and using the unitary invariance of the distance (^), we have 

d{BpkB\ B^+^p^[B^f+^) = 0(2t-^V^) . (59) 

Substituting A: + 1 for /c in Eq. ( |58| ) , we get 

d(p,+i, SpfcSt) = 0(2^('=+i)-^ V^) . (60) 

Using the triangle inequality for the distance measure (^4|), it follows from (|5^) and (^) that 



d(p,+i,B^'+Vo[^^]'+') = 0(2l-^V^) + 0(2i('=+i)-^yf) 

= 0(23('=+y-^^) . (61) 



By induction, this completes the proof of (|55| ) for any k < /cmax- On the other hand 

d(pfc,5Vo[i?^]') = ^Trp|+TVp2_2Tr(pfcSVo[i?t]fc) , (62) 
hence using Eq. ( |55|) it follows that 



/21-'- - 2Tr (pfci?Vo[i?t]'=) = 0(22-- V^) , (63) 
and finally 

IV {p,B''pom') = 2-'[l - 0(-^)] , (64) 



which is equivalent to (4f:) as required. 
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4.3 Atypical initial states 

In this section, we show that the state \0^y^.y'^), where 0^' is a string of r zeros, is an atypical state 
in the sense of the discussion at the end of Sec. ^, i.e., we show that the state \0'^y^.y'^) satisfies 
Eq. (p9|). A direct calculation gives 

1 

Tr[P;^''B\0^y\y^){0^y\y^\B^] = ^ ^ |Ci^*(oVz/', 

£c|=r-l 9=0 

8 '"^'^ 1 + 0(4--") 

- ^ 2^ (4,; + 1)2 • ^^^^ 
11=0 

Substituting t = 2v, we have 



Using Eq. ( p2[ ) and the series representation of Catalan's constant G ~ 0.915965 |26|, 
it follows that 

Tr [P;'i^|0Vy)(0'^y^y'l^t] = ^^^^+0(4'^-")+0(2--) 

~ 0.871 + 0(4'-") + 0(2-'') . (68) 

Since one can treat n — r and r as independent variables, this expression can be made smaller than 
0.872 by choosing n — r and r large enough. For the initial state {O'^y^.y'^), the asymptotic relation 
(|28|) is thus violated. 
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